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Abstract 

We report the results on Big Bang Nucleosynthesis (BBN) based on an updated 
code, with accuracy of the order of 0.1% on 4 He abundance, compared with the 
predictions of other recent similar analysis. We discuss the compatibility of the 
theoretical results, for vanishing neutrino chemical potentials, with the observational 
data. Bounds on the number of relativistic neutrinos and baryon abundance are 
obtained by a likelihood analysis. We also analyze the effect of large neutrino 
chemical potentials on primordial nucleosynthesis, motivated by the recent results 
on the Cosmic Microwave Background Radiation spectrum. The BBN exclusion 
plots for electron neutrino chemical potential and the effective number of relativistic 
neutrinos are reported. We find that the standard BBN seems to be only marginally 
in agreement with the recent BOOMERANG and MAXIMA-1 results, while the 
agreement is much better for degenerate BBN scenarios for large effective number 
of neutrinos, N v ~ 10. 



PACS number (s): 98.80. Cq; 98. 80. Ft 
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1 Introduction 



The synthesis of light nuclei in the early universe, the Big Bang Nucleosynthesis (BBN), 
represents one of the most striking evidence in favour of standard cosmology, and since 
its proposal ]1] , it has been extensively used as one of the best laboratories where to test 
cosmological models and/or elementary particle physics. The appealing feature of BBN 
is that, in its standard version, it relies on quite solid theoretical grounds, which makes 
the predictions for D, 3 He, A He and 7 Li abundances, whose primordial values are only 
partially modified by the subsequent stellar activity, quite robust. 

Recently, the experimental accuracy in measurements of the light primordial nuclide 
abundances, mainly 4 He, has been highly improved, reaching a precision of the order 
of 1%. Similar improvements have been also obtained for both Deuterium (D) and 7 Li 
relative abundances, Yd = D / H and Yiu = 7 Li/H, but, unfortunately, the refinement of 
the experimental techniques does not yet correspond to a clear picture of the primordial 
nuclide densities. This is mainly due to an uncomplete understanding of systematic 
errors. In particular, measuring the primordial A He mass fraction, Y p , from regression to 
zero metallicity in Blue Compact Galaxies, two independent surveys obtained two results, 
a low value 0, 

Y p [l) = 0.234±0.003 , (1.1) 

and a sensibly larger one [[J , 

Y { p h) = 0.244±0.002 , (1.2) 

which are compatible at 2a level only. 

As in a recent analysis [|J, we here adopt a more conservative value, with a larger 
error (hereafter we always use la errors) 

Y p = 0.238±0.005 . (1.3) 

A similar dichotomy holds in D measurements as well, where the study of distant Quasars 
Absorption line Systems (QAS), is thought to represent a reliable way to estimate the pri- 
mordial Deuterium. In this case, observations in different QAS leads to the incompatible 
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results [§, § 1 

Y% ] = (3.4±0.3) 1(T 5 , (1.4) 
Y#° = (2.0±0.5) 1(T 4 . (1.5) 

Finally, a reliable estimate for 7 Li primordial abundance is provided by the Spite plateau, 
observed in the halo of POP II stars [)7|, [8] . The observations give the primordial abundance 
& 

Yr Li = (1.73±0.21) 1(T 10 . (1.6) 

From the theoretical point of view, the BBN predictions are obtained by numerically 
solving a set of coupled Boltzmann equations, which trace the abundances of the different 



nuclides in the framework of standard Big Bang cosmology fllP] . The collisional integrals 



of the above equations contain all n <-> p weak reaction rates and a large nuclear reaction 



network 1 1 



The increasing precision in measuring the primordial abundance has recently pushed 
the theoretical community to make an effort to develop new generation BBN codes [[12]. [L3|] . 
with a comparable level of accuracy. In a recent series of papers, in the framework 
of the standard cosmological model, the present authors |TJJ and other groups [|I2[] have 
performed a comprehensive and accurate analysis of all the physical effects which influence 
4 He mass fraction up to 0.1%. 

Since almost all neutrons present at the onset of nucleosynthesis are fixed into 4 He, 
its abundance is mainly function of the neutron versus proton abundances, at the time of 
n <-> p weak interactions freeze out, which takes place for T ~ 1 MeV. To improve the 
accuracy on the 4 He prediction it is demanding to reach an accuracy level of the order of 
1% in the estimate of the n <-» p rates, well beyond the simple Born approximation. This 
has been performed by considering a number of additional contributions, which, ordered 
according to their relative weight, are the following: 



i) electromagnetic radiative and Coulomb corrections [JL5|, [16 



ii) finite nucleon mass corrections 17, 14 



iii) thermal radiative effects induced by the presence of a surrounding plasma of e ± and 



7 mu 



iv) corrections to the equation of state of the e^, 7 plasma due to thermal mass renor- 
malization jLj, 0; 

v) the residual neutrino coupling to the plasma during e + e~ annihilation, affecting the 
neutrino to photon temperature ratio p(| 0, p~3f| . 

Unfortunately, a similar systematic analysis of all corrections up to the desired level 
of accuracy cannot be performed for the other input parameters of the theory. The 
theoretical estimates do depend in fact on the values of the neutron lifetime, as well as 
on several nuclear reaction cross sections, which are poorly known in the energy range 
relevant for BBN. This introduces a certain level of uncertainty on the light element 
abundances. This aspect has been studied in two different approaches, either using Monte 
Carlo methods [^1] to sample the error distributions of the relevant reaction cross sections, 
or, alternatively, using a linear error propagation p2| , with comparable results. 

In this paper, we further refine our previous predictions by using a new version of our 
numerical code, where the full dependence of the BBN equations on the electron chemical 
potential is accurately implemented. Furthermore, for the standard scenario (vanishing 
neutrino chemical potentials) we perform an accurate likelihood analysis in the space of 
the two free parameters of the model, the effective number of neutrinos, N u , and the 
baryon to photon ratio rj. 

While the value of electron chemical potential is bounded, by neutrality, by the value 
of rj, this is not the case for neutrino-antineutrino asymmetries which, in principle, can 
be quite large. The influence of the neutrino chemical potentials on BBN predictions 
(degenerate BBN) has been considered in the past |[23|| . One relevant aspect of this 
analysis, which is worth stressing, is that degenerate BBN allows for a better agreement 
with observations at values of 77 larger than, say, 6-10~ 10 , while standard BBN prefers 
smaller values, rj ~ (2^6)-10~ 10 . We will discuss this in detail in the paper. This feature is 
particularly relevant in view of the recent analysis of the BOOMERANG and MAXIMA-1 
results on the acoustic peak of the Cosmic Microwave Background Radiation (CMBR), 
which seems to favour a value for the baryonic asymmetry 77 of the order of 10 -9 , larger, as 
we said, than what expected in the framework of standard BBN. This discrepancy may be 
looked as a signal in favour of a large neutrino-antineutrino asymmetry. It seems therefore 
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demanding to re-analyze the degenerate scenario, making profit of the now available more 
precise BBN codes. We have performed this study, and the main result is that degenerate 
BBN and CMBR data seems to be compatible for large values of the effective number of 
neutrinos, N v > 10 and rj < 10~ 9 . 

The paper is organized as follows. In section 2 we briefly describe the BBN set of 
differential equations, recast in a suitable form for a numerical solution. The light element 
abundances obtained from our code, for standard BBN, are then presented in section 3 as 
a functions of r n , N u and 77. Section 4 is devoted to degenerate BBN and to the bounds 
on neutrino chemical potentials coming from nucleosynthesis and CMBR data. Finally, 
in section 4, we give our conclusions. 

2 The BBN set of equations 

Consider N nuc species of nuclides, whose number densities, n^, are normalized with respect 
to the total number density of baryons, n^l] 

Ti ' 

X i = — i = l,..,N nuc . (2.1) 

n B 

Alternatively, we will also make use in the following of the notation: 

X\ = X n , X2 = X p = Xh , X3 = Xd , 
X 5 = Xi He , X 6 = X4 He , X 8 = Xi Li . (2.2) 

The set of differential equations ruling primordial nucleosynthesis is given by 0, [12 



1 = "^m* ■ <-> 

— = -3H , (2.4) 

n B 

p T = -3H(p T + p T ) , (2.5) 
_ / x Nl X Nk X Nl X Nj \ 

Li^Ae) = ^ J2Z jXj , (2.7) 

i 



1 We will use the same notations of our previous paper [O, which we refer to for further details. 
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where and pr denote the total energy density and pressure, respectively, 

Pt = p 7 + Pe + Pu + Pb = Pnb + Pb , (2.8) 
Pt = P 7 + Pe + Pu + Pb , (2.9) 

i,j,k,l = (1, .., N nuc ), Zi is the charge number of the i— th nuclide, and the function 
L(z,y) is defined as 

X roo / 1 1 \ 

L(z,y) = — dxxVx 2 -z 2 . (2.10) 

y y> ir 2 Jz \e x ~y + i e x +y + \) y ' 



Eq.( |2.3| ) is the definition of the Hubble parameter, H, whereas Eq.s ( |2.4j ) and ( j2.5| ) state 
the total baryon number and entropy conservation in the comoving volume, respectively. 
The set of N nuc Boltzmann equations Q2.6Q describe the density evolution of each nuclide 
specie, and finally Eq. (|2.7| ) states the universe charge neutrality in terms of the electron 
chemical potential, <ft e = /x e /T, with T the temperature of e , 7 plasma. Note that the 
neutrino energy density and pressure are included in Eq. (|2.5|) only for T > Tjj^ 



In a previous analysis, [13], we neglected the contribution of e in the BBN equations, 
since its effect results to be very small. In this way we obtained a substantial simplification 
of the set of equations ( |2.3|) - (|2.71) , since the unknown functions can be reduced in this case 
to the N nuc +1 (h = ub/T 3 , Xj). We do here release this assumption and report the results 
of an improved code where we take into complete account the electron chemical potential 
evolution. We will see, however, that all changes on the final abundances are of minor 
impact. In order to obtain the new equations, we note that it is more convenient to follow 
the evolution of the N nuc + 1 unknown functions (0 e , Xj) in terms of the dimensionless 
variable z = m e /T, and to use Eq. (|2.7|) to get % as a function of <p e . The new set of 
differential equations may be cast in the form 

d(j) e 1 L E F + {z L z - 3 L) G 



dZ Z LEtfc-L+G 
dX< L, e F + (zL z -3L)lf e 

dZ ^ '* LE lt- L *< G 



2.11 



(2.12) 



2 We assume that all neutrinos decouple at the same temperature, To — 2.3 MeV 
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where the functions E, F and G are given by 

£ ZiTi 



E{z, <j) e , Xj) 



F(z, e , X 3 ) 



G(z, e , X 3 ) 



3H 



H Z 3 X 3 



3 A 5p e 

2 02 

3if + p) + S E (am* + ^ ) r, 



(2.13) 

(2.14) 
(2.15) 



-<3 3 3 i 

and H = m e H, n B = ml n B , Tj = m e fj, p T = T 4 p T , p T = T 4 p T . Note that by using 
Eq. ( J2.7|) and the previous definitions, it is possible to express ps and p B as functions of 

z, 4> e , and Xj only 



PB 



PB 



z Liz, 



(2.16) 
(2.17) 



With AMj and M u we denote the i-th nuclide mass excess and the atomic mass unit, 
respectively, normalized to m e . In Appendix A we report the partial derivative of L with 
respect to z and e , denoted with L z and L^ e , and the quantities p e , 5p e /5z and 5p e /5<fi e 
in a form which is suitable for a BBN code implementation. 

Eq.s (gHD - ([2l2D are solved by imposing the following initial conditions at Z{ n = 
m e /(10MeV): 



<\>e(z in ) = e ° 



(2.18) 



Xi(Zin) = f C(3)^ 



(exp{g^ n } + 1) 1 , X p (z in ) = (exp{-q z in } + 1) 1 , (2.19) 

Aj-l 3 

/ p n 



A} 



MmZ, 



x exp 



with z = 3, ..,N nuc 



[2.20) 



In the previous equations g = (M n — M p )/m e , and the quantities A { and denote the 
atomic number and the binding energy of the i— th nuclide normalized to electron mass, 
respectively. Finally r\ is, as usual, the baryon to photon number density ratio, and (p e ° 
the solution of the implicit equation 



LyZi n j 



11 2C(3) 

4 7T 2 



7/ £^Xj 



(2.21) 



The method of resolution of the BBN equations ( |2.11| )-( p.l2j ) is the same applied in [ |D| . 



It is the Backward Differentiation Formulas with Newton's method, implemented in a 



NAG routine with adaptive step-size (see |T3[ for more details). We used the reduced 



network of nuclear reactions, made of 25 reactions involving 9 nuclides, since the use of 
the complete network affects the abundances for no more than 0.01 %. 

3 Primordial abundances for standard BBN 

The new BBN code has been used to produce the primordial abundances, in the standard 
scenario, for different values of the input parameters, namely the neutron lifetime r n , the 
effective number of neutrinos N v , defined as 

with p v the total neutrino energy densities, and the final baryon to photon number density 
ratio n. 

We consider the abundances relative to hydrogen for D, 3 He, and 7 Li, 

Y D = ^ , Y He = ^ , Yr u = ^ . (3.2) 

A. H A# A. H 

In the case of 4 He, the quantity usually defined as mass fraction, 

Ai u„ Xi Up , 

Y P = " A X ' ^ 
L,j A i A i 



is rather the baryon number fraction, since A± He = 4. Note that expression (|3.3|) does not 
correspond to the true mass fraction, obviously defined as 

y m = Mi He Xi He 

" Ej M, x 3 ■ 1 • > 

The difference between Y p and Y™ are of the order of 1% and thus relevant for an accurate 
analysis like the one presented here. Since it is customary to express the experimental 
value in terms of Y p , see (|1.3| ), we will consider this quantity for a comparison with 
experimental data. 

Table |1] shows the results obtained for N„ = 3, r n = 886.7 s, and rj = 5-10~ 10 , compared 
with our previous results [ 13| . As one can see, the inclusion of the complete evolution of 
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Y D 


Ys He 




"yrn 
V 




present 
analysis 


0.3609-10" 4 


0.1166-10" 4 


0.2461 


0.2448 


0.2879-10~ 9 


previous 
analysis 


0.3638-10- 4 


0.1175-10- 4 


0.2460 


0.2447 


0.2814-10- 9 















Table 1: The predictions on light element abundances obtained with the new BBN code, 



compared with our previous results |L3 



e modifies the final A He abundances for less than 0.1% and for few percent the other 
nuclides. 

In Figures [1], [| |^, [7] we show the theoretical predictions for the abundances (|3.2|)-(|3~3|). 
compared with the experimental values. Figures |2|, f|, |6], |8| show the comparison between 
our results and the ones of Ref.s [0, [12]. The relative differences among our predictions 



and the results of Ref.s 12] for 4 He, in the relevant range for 7710 = 10 10 ^, are less 
than 0.25%, and thus probably due to different ways of taking into account subdominant 
effects (thermal radiative corrections). The differences are not much larger for Deuterium, 
but reach few percents for 3 He and 7 Li, which however have very large theoretical errors 
due to the uncertainties on nuclear reaction rates. Note that for D, 3 He and 7 Li only the 
results of |22| are available. 

We have performed a fit of the previous abundances as functions of N v , r n , and 
x = log 10 r] 10 , with an accuracy which is better than 1% in the all ranges 0.5 < N v < 6, 
882.9 s < r n < 890.5 s, < x < 1. In particular, for N v = 3 and 4 He, our fit is accurate 
at 0.01%. The fitting functions have been chosen as 



ki ■ Y{ (iVy, T n , x) 



j^a^+lj^bjxA (Tn-O+fEci^) (iV,-3) + 

3=0 \j=0 J \j=0 J 

E d J xj j (N„ - 3) 2 ] exp I e 3 x 3 J , (3.5) 



with all coefficients given in Table [| 

By using these expressions ( |3.5| ), as the theoretical predictions for the light element 
abundances, and their experimental measurements (|1.3|) , ( |i~4]) , ( p..5|) , (|1.6| ) it is possible 



1 XT 


J 


a i 




C 3 


dj 




rm — ~ 1 

W-D 





0.48212 


- 


0.076014 


0.0015595 


- 




1 


0.12998 


- 


0.011704 


0.0015438 


-4.0240 




2 


0.24279 


- 


-0.056835 


-0.010037 


-1.1394 




3 


-0.91776 


- 


0.0045225 


-0.014180 


4.7712 




4 


2.2660 


- 


0.27257 


0.090586 


-9.7616 




5 


-1.1308 


- 


-0.20773 


-0.15266 


9.5377 




6 


-1.4315 


- 


-0.096193 


0.11416 


-3.5489 




7 


1.0751 


- 


0.10664 


-0.031562 


- 


10 5 - 3 He 





3.3201 


- 


0.18342 


-0.0067466 


- 




1 


-8.8146 


- 


-0.60698 


0.013121 


0.93721 




2 


7.3100 


- 


1.0607 


-0.042731 


2.1792 




3 


14.543 


- 


-0.52115 


0.26524 


-5.6965 




4 


-51.298 


- 


-1.9291 


-0.72719 


5.8729 




5 


76.092 


- 


4.5046 


0.94031 


-8.2579 




6 


-59.309 


- 


-3.8725 


-0.58402 


3.8572 




7 


20.591 


- 


1.2760 


0.13864 


- 


10- 4 He 





2.2289 


0.0020479 


0.13021 


-0.0096485 


- 




1 


-0.052824 


-0.00064726 


0.0075762 


0.0018598 


0.27365 




2 


0.59320 


0.0033813 


-0.013789 


-0.0030398 


-0.58807 




3 


0.49399 


-0.014272 


0.082061 


-0.011842 


0.22190 




4 


1.9579 


0.041398 


0.10938 


0.027656 


-1.3763 




5 


-0.89742 


-0.057323 


-0.12581 


-0.051066 


0.93037 




6 


0.34269 


0.040258 


0.096350 


0.039127 


-0.17224 




7 


0.47527 


-0.010753 


-0.0018617 


-0.013441 


- 


10- A He m 





2.2184 


0.0020421 


0.13003 


-0.0095937 


- 




1 


-0.44715 


-0.0011625 


-0.016604 


0.0040428 


0.44707 




2 


0.94939 


0.0065040 


0.010744 


-0.013907 


-0.69814 




3 


6.8365 


-0.026617 


0.43785 


0.018033 


-2.8141 




4 


-1.1513 


0.094984 


0.0086552 


-0.12409 


-0.80120 




5 


2.0116 


-0.14739 


-0.17116 


0.19005 


3.0861 




6 


3.1081 


0.11736 


0.48418 


-0.16606 


-1.1642 




7 


4.0381 


-0.031739 


0.088568 


0.031865 


- 


10 9 - 7 Lz 





0.52920 


- 


0.15387 


0.011486 


- 




1 


-1.5617 


- 


-0.42097 


-0.0075472 


-2.0190 




o 
Z 


z.uuuz 




-U.ZoUOo 


-U.Zlo ( o 






3 


0.58298 




2.9150 


0.94226 


9.8088 




4 


-8.3304 




-6.7141 


-1.7403 


-9.7328 




5 


22.627 




6.2727 


1.8244 


3.1404 




6 


-21.053 




-1.9273 


-1.1338 


0.15166 




7 


6.1684 




-0.090045 


0.32046 





Table 2: Values of the coefficients of Eq. ( |3.5|) for light element abundances. 
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to test the compatibility of standard BBN scenario in the N u —T]iq plane. To this end, we 
define a total likelihood function as 



£{Nv, Vw) = l d{N u , 7710) L4 He (N u , rjio) L7 Li (N v , rj 10 ) , (3.6) 

where the likelihood function for each abundance, assuming Gaussian distribution for the 
errors, is given by the overlap 

Li(N v , 7710) = 

1 r dY c::p f ( y - ^q)) 2 l c::p f ( Y - y n 2 } ( o 7) 

In order to evaluate Li(N v , rjio) we need the theoretical uncertainties aj h (N u , r] W ). In 
Ref. |22| a new, alternative approach to the standard Monte Carlo technique [^l|] has 
been proposed. This method is based on linear error propagation for the estimate of 
theoretical uncertainties on primordial abundances due to the poor knowledge of nuclear 
reaction rates. The light element abundances, Yj, and the logarithmic derivatives of Y{ 
with respect to the nuclear rates, A^, are given as polynomial fits, while the variation of 
the abundances for a change 5Rk in the rate Rk is obtained as 

5Y t = Y l Y. A,. fc ^ • (3.8) 

k R k 

Taking into account the error correlation, the error matrix results 

°ij = Yi Y i H XikX i k i~R^) ' ( 3 - 9 ) 

where AR^ are the la uncertainties. In particular the theoretical la uncertainties are 
given by the square root of the diagonal elements, 

th _ f~2 



We used the Fortran code provided by the authors [p6] to calculate the theoretical uncer- 
tainties af l {N v , 7710), which we used in Eq. ( |3.7|) . 

The total likelihood function (|3.6| ) and the corresponding contour plots for 50%, 68% 
and 95% CL, for low and high D are shown in Figures |9], |H] and [11], [12|, respectively. 
As already clear from Figure [I], the two different experimental estimates of D single out 

11 



different regions for 7y 10 . From the 95% CL contour of Figure [10], for low D we have 



4.0 < rjio < 5.7, which is in fair agreement with the similar results of Ref. Q], whereas 
for high D from Figure |12| we get 1.4 < i] w < 3.7. 

As far as N u is concerned, for low D and for 95% CL one has 1.7 < N v < 3.3, whereas 
for high D one gets 2.3 < N v < 4.4. In both cases we obtain comparable ranges for N v 
with respect to Ref. Q]. The total likelihood function is peaked around the points shown 
as the crosses in Figures [II] and [L2|, which correspond to N u = 2.44 and r)io = 4.69 for 



low D, and to N u = 3.29 and rjio = 1.81 for high D. The position of the two maxima 
is easily understood. The Deuterium abundance is a decreasing function of 7710- Since 
lowering 7] w results in a smaller 4 He mass fraction Y p , it is necessary to compensate this 
effect by increasing the universe expansion rate via a larger N v . This in fact leads to a 
larger value for the freeze-out temperature for nucleon weak interactions and thus gives 
a larger amount of the initial neutron to proton density ratio. The single contributions 
( |3.7| ) to the total likelihood ( |3.6| ) can be easily recognized by looking at Figures [13| - p~6 



where we show the D, A He and 7 Li likelihood functions, for both high and low Deuterium 
results, for the two preferred values N v = 2.44 and N v = 3.29. For the low D case, a 
better overlap between the maxima of the single likelihoods of D and 4 i7e is realized for 
N v = 2.44. The opposite situation occurs for high D where to N„ = 3.29 corresponds a 
better overlap of the single likelihoods. 

4 BBN predictions for degenerate neutrinos 

The assumption of vanishing neutrino chemical potentials can be only justified by the 
sake of simplicity or by the theoretical prejudice that their order of magnitude should 
be set by the ones of the corresponding charged leptonic partner. However, physical 
scenarios in which large lepton asymmetries are produced, which do not lead to a large 



baryon asymmetry, have been proposed in literature [g7], gg, ^9], |3(|. They are based 



on the Affleck-Dine mechanism [^7| or on active-sterile neutrino oscillations |3T|, In 



particular, the expected asymmetry may be different for each neutrino family. For these 
reasons, at least in principle, it is worth-while considering the primordial nucleosynthesis 
in presence of large neutrino asymmetries, i.e. for non-vanishing neutrino chemical po- 
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tentials. It seems to us that this topic, which has been extensively studied in past [ f23|| , 
receives a renewed interest in view of the recent BOOMERANG and MAXIMA- 1 results 



24j| on the acoustic peaks of the CMBR, which suggests a larger value of the baryonic 
matter contribution to the total energy density Q. In their analysis of the data they find 
a baryon density Vt^h 2 ~ 0.02 -r- 0.03, or rj ~ (6 4- 10) 10~ 10 , which, as is clear from the 
considerations of the previous section, is incompatible at 95% CL with the standard BBN 
result (see Figures |10| and [12]). Actually even larger values for VL^h 2 are obtained if no 



constraints are imposed in the likelihood analysis (see [plf). It is the aim of this section 
to discuss whether a finite neutrino chemical potential may reconcile BBN theoretical 
predictions, the observed nuclide abundances and a higher value for rj. In particular we 
report the results of a new analysis of the degenerate BBN scenario we have performed 
with our code. 

The effect of neutrino chemical potentials on BBN predictions is twofold. Due to 
the definition of N v (|3.1|), non- vanishing £ a = fi Ua /T u , with a denoting the neutrino 
specie, change its value from the non-degenerate case (£„ = 0). In fact for three massless 
neutrinos with degeneracy parameter £ Q , N v becomes 



/v — "\ -\- y, 

ly u i /L - , a=e,[i,T 



(4.11) 



implying a larger expansion rate of the universe with respect to the non-degenerate sce- 
nario. In this case, nucleons freeze out at a larger temperature, with a higher value for the 
neutron to proton density ratio, which implies a larger value of Y p . This effect does not 
depend on the particular neutrino chemical potentials but rather on the whole neutrino- 
antineutrino asymmetry, via the sum in the r.h.s. of (|4.11|) . In addition, we also note that 
the neutrino decoupling temperature and the ratio T v /T, entering in the BBN equations, 
are affected by a change of N u as well. However, it has been checked that this effect 
is quite negligible on the predictions for the element abundances. 

Electron neutrinos entering in the n <-> p processes, can modify the corresponding 
rates if their distribution has a non- vanishing £ e . In particular, a positive value for £ e 
means a larger number of v e with respect to v e and thus enhances n — > p processes with 
respect to the inverse processes. This, of course, reduces the number of neutrons available 
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at the onset of BBN. Moreover, since the initial condition for BBN at T ~ 10 MeV is fixed 
by the nuclear thermal equilibrium, which is also kept by reactions like v e + n — > p + e~, 
the chemical equilibrium fixes — fi n )/T ~ £ e . This implies an initial n/p ratio which 
is lowered by the factor exp(— £ e ), which again reduces the number of neutrons available 
at the onset of BBN. 

These effects strongly influence the nuclide production, so with no and v T degener- 



acy, the value of £ e is strongly constrained. In this case, the authors of |23[ found the limits 
—0.06 < £ e < 0.14. For a fully degenerate BBN, at least for Y p , the effect of a positive 
£ e can be compensated by the contribution to N v coming from £ MjT . In |23| the neutrino 
degeneracy parameters result to be in the ranges —0.06 < £ e < 1.1 and |£ MiT | < 6.9 for 
7^ 0, ir = or viceversa, and \^, T \ < 5.6 for £ M = £ T 7^ 0. 

We do consider in our analysis as input parameters £ e , N v and 7710- Likelihood analysis 
of compatibility between theoretical predictions and experimental data yields contour 
levels which are surfaces in the three-dimensional space of these parameters. Because of 
the partial cancellation of the effects due to £ e and £ MiT we have discussed, if we define, 
in analogy to the non-degenerate total likelihood function, £(£ e , N„, r] w ), it is 

reasonable to expect that it may sensibly differ from zero in a quite wide parameter region, 
so some bound should be chosen to the possible range of variation for the parameters. 
We have chosen to constrain N v to be smaller than 13. This bound has been obtained, at 



2a level, in [34|, by a likelihood analysis of the BOOMERANG data, as function of N u , 
maximizing for each N v the likelihood function over all other parameters, including r/. As 
we will see our conclusion on the degenerate BBN scenario versus BOOMERANG and 



MAXIMA- 1 data is completely different than what has been argued in [p4| . Nevertheless 



we think that the upper bound on N v is quite robust. We also consider £ e and r] w in the 
wide range —1 4- 1 and 1 30, respectively. 

By using the results of our BBN code for the degenerate neutrino case and performing 
a study similar to the one presented for the standard scenario, we obtained likelihood 
functions for both low and high D experimental values. We first report the maxima for 
these functions in the low D case 

£ e = 0.06 , N u = 3.43 , tii = 5 , (4.12) 
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and in the high D scenario 

£ e = 0.35 , jV„ = 13 , ^ = 4.20 , (4.13) 

We notice that in the low Deuterium case the solution prefers an almost non degenerate 
scenario, with values for 77 compatible with the one obtained in the previous section and 
a slightly larger result for N v . 

More interesting is to follows the allowed ranges for £ e and N u as functions of 77. In 
Figures [T7| and [18] we show £(£ e , A^, 7710), evaluated for 7710 taking the values in ( |4.12| ) 



and (|4.13| ). Increasing 7710 in the considered range, the likelihood functions move towards 
higher values of £ e and N v , along, approximatively, a linear path. This is due to the fact 
that increasing £ e , which results in a lower value for the n/p ratio, and then for Y p , must 
be compensated by a faster expansion produced by a higher value for iVj,. This behaviour 
is highlighted in the 95% exclusion plots for the £ e and N v parameters for different values 
of 7710, which are reported in Figures [19] and |20] with the bound N v < 13. These two plots 
summarize our analysis of the degenerate BBN scenario, providing the combined bounds 
on £ e and N v . For low D the allowed range for i] w is 3.3 4- 9.9, while for high D we have 
1.1-r- 5.8. In both cases the degenerate BBN scenario, at 2a level is compatible with 
observational results on nuclide abundances even for quite large values for both N„ and 
£ e . As expected, however, these two parameters are strongly correlated. 

From our analysis we see that a large value for 77 ~ 10 -9 is compatible, at 2a level, 
with degenerate BBN for quite large N v only, N v > 10, and large positive £ e > 0.3. The 
BOOMERANG and MAXIMA-1 results seem to indicate a value for 77 ~ 10~ 9 , but at the 



same time suggest N u < 13, again at 2a level [34]. We may conclude that the degenerate 



BBN scenario and BOOMERANG and MAXIMA-1 results are compatible only in the 



very top area in the exclusion plot of Figure (|T9|) for the low D scenario. In the high D 
case, the upper limit on N v is already reached for 77 ~ 6.5- 10~ 10 . As a conclusion we can 
say that, if the result 77 ~ 10~ 9 will be confirmed, the degenerate BBN scenario is not 
consistent with CMBR results at 2a level, while there is agreement for the low D case. 
For slightly lower values for 77 ~ (0.5 -j- 0.6) 10 -10 , which corresponds to the central value 
measured by MAXIMA-1 experiment, the agreement at 2a level definitely improves, and 
it would represent a signal in favour of a large neutrino-antineutrino asymmetry in the 
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universe, 0.2 < £ e < 0.5 and N v > 10 P5] . We note that a similar result has been obtained 
in 



5 Conclusions 

In this paper we have reported the results of a likelihood analysis of Big Bang Nucleosyn- 
thesis theoretical predictions versus the experimental data, both for vanishing neutrino 
chemical potentials (standard BBN) and in presence of neutrino-antineutrino asymmetries 
(degenerate BBN). The theoretical estimates have been obtained by using a new updated 
code we have developed in the recent years to increase the accuracy at the 1% level. 

In the standard scenario, BBN predictions are in good agreement with experimental 
data for both low and high D results. In the first case, since Yd is a decreasing function 
of baryonic asymmetry 77, the low Deuterium abundance fixes the maximum of total 
likelihood to larger values of 77, and to compensate the effect on Y p yields smaller values 
of N v . In this way we have obtained as preferred values N v = 2.44 and 7710 = 4.69. 
The situation is reversed for high D where the maximum lies at N v = 3.29 and r] w = 
1.81. From the plots of 95% CL for both the above likelihood functions one gets the 
corresponding compatibility regions 

low D 1.7 < N v < 3.3 4.0 < r) W < 5.7 0.015 < Vt B h 2 < 0.021 

. (5.14) 

high D 2.3 < N v < 4.4 1.4 < r) 10 < 3.7 0.005 < Q B h 2 < 0.014 

Both ranges on N u are comparable with the results of Ref . [|J . For the baryon asymmetry 
T], the result for low D fairly coincides with whereas the distortion of the contour for 
95% CL in case of high D makes our upper limit on r]io larger than the value of ||J]. 
This effect is mainly due to the non-trivial dependence on the parameters N v — r] w of the 
nuclear theoretical uncertainties. 

We have also analyzed the degenerate BBN scenario, motivated by the BOOMERANG 
and MAXIMA- 1 results |24| on the CMBR spectrum. For both high and low Deuterium 



measurements we have obtained a 2a exclusion plots for the £ e and N u parameters, and 
we have shown that the theoretical estimates are in good agreement with the experimen- 
tal measurements of light nuclei abundances even for large values of these parameters, 
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provided they lie in the regions shown in Fig.s ( |T9"D and (|20|). In particular large values 
for 1] implies high values for the neutrino degeneracy parameters. 

A first analysis of the recent CMBR data, f24|, suggests quite large values for r\ ~ 1CT 9 . 
From our results we see that this baryon to photon density is actually too high to be in 
good agreement with degenerate BBN, the compatibility being only at 2a level, with the 
low Deuterium scenario. This conclusion holds if we constraint the effective number of 



neutrinos to be bounded by N u < 13, which, as pointed out in |34| , is again suggested 
by BOOMERANG data. The agreement improves for slightly smaller values for 77, as the 
central value obtained by MAXIMA-1 results, and would be an evidence in favour of large 
neutrino degeneracy, 0.2 < £ e < 0.5, N v > 10, an intriguing feature of the neutrino cosmic 
background. New experimental results, as the one expected from the MAP and PLANCK 
experiments on CMBR, as well as new analysis of the BOOMERANG and MAXIMA-1 
data, will be crucial in clarifying this issue. 
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A Series expansion for integral quantities 

The partial derivatives of L with respect to z and <p e , L z and L^ e , can be expressed as 
the following series^ 

L z (z, e ) = ^ \ — £ (-!) n+1 [Ki(n z R ) + K 3 (n z R )] sinh(n0 e ) , (A.l) 

n=l J 

2 4 °° 



dz { Zr V 71 J 

L^{z, e ) = ^ J2 (-!) n+1 K 2 {nz R ) cosh(n0 e ) , (A.2) 



* n=l 



where the Ki are the modified Bessel functions of the second kind and in z R (z) = m^/T 
we consider the renormalized electron mass, 



(T) ~ m e 



71 




l + - a 






\m e J 



(A.3) 



3 We truncate all the series at n — 7 while in the standard code jl0| the truncation is at n = 5. 
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The series expansion for the remaining electron quantities are 



3 oo (_-\\n+l 

p e (z,<j) e ) = — [SK^nzRl+KtinzR)} cosh(n</» e ) , (A.4) 

^ ^ n=l n 



^,0e) = ^{^ + ^E (-l) n+1 



cosh(n0 e ) \ , (A. 5) 



5p e z; 



3 oo 



.2 



) fey 



2 7T 2 

z 71 n=l 



£ (-l) n+1 [3^3^^) +^(71^)] sinh(n0 e ) . (A.6) 
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Figure 1: The Deuterium abundance Yd (|3.2|), for N v = 3, versus 7710 = 10 77. The hori- 
zontal bands correspond to the experimental determinations of Ref.s [||, ||. The solid lines 
bound the theoretical predictions at la t h- The dashed line, on this scale indistinguishable 
from the central solid line, is the result from p2|. No result is given for D in [|12|. The 



same notation has been applied in the following Figures. 
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Figure 2: Relative difference (in percent) for Deuterium abundance between the results 
of the present work and the ones from El in the case N u = 3. 
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Figure 4: Relative difference (in percent) for Y-$He between the results of present work and 
the ones from [22] in the case N v = 3. 
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Figure 7: The lithium-7 abundance Y7 Li ( |3.2| ), for N„ = 3, versus r] W . The dashed line, on 



this scale indistinguishable from the central solid line, is the result from [22]. No result is 



given for 7 Li in [12 
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Figure 8: Relative difference (in percent) for Y7 Li between the results of present work and 



the ones from 213] in the case N v = 3. 
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Figure 10: Contour plots of the total likelihood function for the low experimental value 
of D. From inner to outer they correspond to 50%, 68% and 95% CL, respectively. The 
cross indicates the maximum of likelihood function, and corresponds to N v = 2.44 and 
Via = 4.69. 
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Figure 12: Contour plots of the total likelihood function for the high experimental value 
of D. From inner to outer they correspond to 50%, 68% and 95% CL, respectively. The 
cross corresponds to N u = 3.29 and r] W = 1.81. 
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Figure 13: Single likelihood functions for the three light abundances in the case of the 
low experimental value of D and N u = 3.29. The single peaked curve stands for Ld, the 
double peaked curve for L7 Li and the broad one for L4 He . The same notation holds for 
the following figures. 
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Figure 15: Single likelihood functions for the three light abundances in the case of the 
high experimental value of D and N u = 3.29. 
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Figure 17: Total likelihood function in the degenerate BBN scenario in the plane £ e -iVj 
in the case of the low experimental value of D for rjio = 5. 
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Figure 18: Total likelihood function in the degenerate BBN scenario, in the case of the 
high experimental value of D for i] 1Q = 4.20. 
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Figure 19: The 95% exclusion plot for the variables £ e — N v for low D experimental value. 
The dark area represents the BBN allowed region, for 7710 in the range 3.3^-9.9. The cross 
corresponds to the maximum of the likelihood function. 
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Figure 20: The 95% exclusion plot for the variables £ e — N v for high D experimental 
value. The dark area represents the BBN allowed region, for rjio in the range 1.1 -j- 5.8. 
The cross corresponds to the maximum of the likelihood function. 
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